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Abstract 

A new method for solving the time-dependent two-center Dirac equation is developed. The time-dependent 
Dirac wave function is represented as a sum of atomic-like Dirac-Sturm orbitals, localized at the ions. The 
atomic orbitals are obtained by solving numerically the finite-difference one-center Dirac and Dirac-Sturm equa- 
tions with the potential which is the sum of the exact reference-nucleus potential and a monopole-approximation 
potential from the other nucleus. An original procedure to calculate the two-center integrals with these orbitals 
is proposed. The approach is tested by calculations of the charge transfer and ionization cross sections for the 
H(ls)-proton collisions at proton energies from 1 keV to 100 keV. The obtained results are compared with related 
experimental and other theoretical data. To investigate the role of the relativistic effects, the charge transfer cross 
sections for the Ne^+(ls)-Ne^°+ (at energies from 0.1 to 10 MeV/u) and U^^+(ls)-U^^+ (at energies from 6 to 
10 MeV/u) colhsions are calculated in both relativistic and nonrelativistic cases. 

PACS numbers: 34.10.+X, 34.50.-s, 34.70.+e 
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I. INTRODUCTION 



Since the pioneering works |lll43|], where the oscillatory behavior of the resonance charge-transfer 
probability for low-energy collisions was predicted, numerous publications have been devoted to 
the theoretical investigations of the charge-transfer, excitations and ionization in the H(ls)-H+ col- 
lisions (see, e.g., review s p -tal). Nonrelativistic two-center finite basis set calculations have been 
carried out in Refs. OKISD. Nonrelativistic three-dimensional lattice methods in the position and 
momentum spaces have been applied for the time-dependent Schrodinger equation in Refs. [|l4l- ll8|l . 
Within the nonrelativistic approach, the probabilities and cross-sections for a homonuclear collision 
for the nuclear charge Z > 1 can be easily obtained by scaling to the H(ls)— H+ 
collision. In the straight- line trajectory approximation, the cross section a(Z,v) scales exactly as 
a{Z,v) = a{l,v/Z)/Z'^ [fl [3], where v is the projectile velocity. This scaling law is not valid, 
however, in the relativistic theory. 

Collisions involving highly char ged ions provide tests of relativistic and quantum electrodynamics 
effects in the scattering theory ll20l - l22ll . The study of such processes can provide also a unique tool 
to probe the quantum electrodynamics (QED) in the supercritical Coulomb field, if the total charge of 
the colliding ions Z = Za + Zb is larger than the critical one Zc = 173 (see, e.g., Refs. [|23l-l26l] 
and references therein). In the presence of such a field the energy of the one-electron 1(t+ state of the 
quasi-molecule can reach the negative-energy Dirac continuum, when the distance R between target 
ion A and projectile ion B becomes equal to the critical value Rc. For the distance R less than Rc 
the ground state level dives into the negative-continuum spectrum. In the U^^+(ls)-U^^+ collision the 



critical radius for the point nucleus case was found to be Rc = 36.8 fm [ |27l l. 

To date various approaches were developed to treat the heavy-ion collisions [|28l1 . In Refs. [129 
the two- and three-dimensional numerical lattice methods were employed to solve the time-dependent 
Dirac equation at high energies. In Refs. OVi - koll . high energy relativistic collisions of heavy ions were 
considered using the basis set approach, in which the time-dependent wave function was expanded in 
terms of the atomic eigenstates of the projectile and the target. For internuclear distances smaller than 
about 1000 fm some effects can also be evaluated within so-called monopole approximation, which 
accounts only for the spherically-symmetric part of the two-center potential [|4ll-l43ll. The atomic pro- 
cesses such as excitation, ionization and charge transfer in relativistic atomic collisions involving heavy 
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and 



lighly-charged projectile ions with energies ranging from 100 MeV/u upward were studied in Refs. 



34, 



44l- l46|l and references therein. 



In the present work, we develop a new method for solving the two-center stationary and time- 
dependent Dirac equations. The wave functions are expanded in terms of the Dirac and Dirac-Sturm 
basis functions, which are central-field 4-component Dirac bispinors centered at the ions. The ra- 
dial parts of these orbitals are obtained by solving numerically the finite-difference radial one-center 
Dirac and Dirac-Sturm equations. In the nonrelativistic calculations of atoms and molecules, so-called 



Coulomb-Sturmian basis set was introduced in Ref. 1471] . The Hartree-Fock calculations of atoms with 



this basis were considered by many authors (see, e.g., Ref. [|48|]). The relativistic Coulomb-Sturmian 
basis was employed in the papers [|49l- l52|l . In the present paper we use a non-Coulomb relativistic 
Sturm basis set, which is obtained by solving numerically the Dirac-Sturm equations with a special 



choice of the weight function, that was proposed in Refs. [|53 



5411 . This allows us to include any 



central-field potential in the radial equations for the large and small components of the basis functions. 
In particular, the Coulomb potential of the other ion can be included in the radial equations within the 
monopole approximation. The basis set constructed in this way is described in detail in section HTBI 
Calculations of two-center integrals with the basis functions obtained require using special tools. In 



the nonrelativistic case, a special symmetrica 



OTOcedure for such calculations, based on the Lowdin 



57h . In section HTCl we generalize this procedure to the 



reexpansion [15 511 . was developed in Refs. ||5 
relativistic case. 

To test the quality of the two-center expansion described above we perform relativistic calcula- 
tions of the ground-state energy of molecular ion and one-electron quasi-molecule Th2'^^^ at the 
"chemical" distance R = 2/Z a.u. and compare the results with high-precision calculations of Refs. 

We also calculate the ground-state energy as a function of the intemuclear distance R and the 
critical radii Rc for a number of one-electron quasi-molecules, including U2*^^. Most calculations of 
the critical distances Rc presented in the literature were performed either for the point-nucleus model 



601, 



611 ] or with a crude estimate of the nuclear-size effect ^^M- We calculate the critical dis- 



tances for both point and extended nucleus models using the same basis set expansion. The obtained 
results and comparison with the calculations by other authors are presented in section HITBI 

The classical Rutherford trajectories n25\\ of the projectile and target ions are obtained by numerical 
solution of the Newton's equations. The Bom-Oppenheimer approximation is used to separate the mo- 
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tion of the electron and the nuclei. The magnetic interaction between the electron and the moving ions 
is neglected, because of low velocity of the projectile with respect to the target. The time-dependent 
Dirac equation for the electron is solved using the two-center basis set expansion. The expansion 
coefficients can be defined employing, e.g., the Crank-Nicholsen propagation scheme [65] or the split- 
operator method [|66n . These methods conserve the norm of the time-dependent wave function at each 
time step, since the Crank-Nicholsen operator and the split-operator are unitary. However, in this work 
we use the direct evolution exponential operator method, which is more stable compared to the others. 
To obtain the matrix of the exponential operator in the finite basis set one has to diagonalize the gener- 
alized Hamiltonian matrix at each step of time. Since our basis set is not too large, the diagonalization 
procedure is not too time consuming. The amplitudes of the charge transfer to different bound states of 
the projectile ion are calculated nonperturbatively by projecting the time-dependent wave function onto 
the moving Dirac orbitals of the projectile. 

In section IIII CI we present the results of the relativistic calculations of the charge-transfer proba- 
bilities and cross sections for the H(ls)-H+, Ne9+(ls)-Nei°+, Xe^^+(ls)-Xe^^+ , and U^^+(ls)-U^'^+ 
low-energy collisions. All the calculations are performed in the laboratory frame S, that is defined to 
be at rest with respect to the initial target position. The H(ls)-H+ collision is considered in section 
IIII C 1 1 Since the relativistic effects in this collision are negligible, the results of our calculations can be 
compared with nonrelativistic data obtained by other authors ( section UlIC II) . The role of the relativis- 
tic effects is investigated in sections HlIC 2llIlIC 31 and llllC 4l where the relativistic and nonrelativistic 
calculations of the charge-transfer probabilities and cross sections are performed for higher- Z ions. 



II. THEORY 



A. Two-center Dirac equation in the finite basis set 



1. Two-center expansion 

Within the Born-Oppenheimer approximation, the motion of the electron is considered as a mo- 
tion in the field of the two nuclei being at given positions (the stationary case) or moving along the 
classical trajectories (the non- stationary case). Let Ra and Rb are the positions of the target {A) and 
projectile {B) nuclei, respectively. The time-dependent \E'(r, t) and stationary il){r) wave functions are 
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the solutions of the time-dependent and stationary Dirac equations, respectively. In the atomic units 
(h = m = e = 1), these equations are given by 

(9\&(r t) 

i — = /iD^(r',t), /iDV^„(r) =e„^„(r). (1) 
Here is the energy of the stationary state and the /id is the two-center Dirac Hamiltonian defined by 

hj) = c{cfp) + /3c^ + VAB{r) , (2) 
where c is the speed of light, a., (3 are the Dirac matrices, and 

VabW = KuciM + Kuci(^b) , rA = r-RA. rj, = r - Rb, (3) 
-Z/r for the point nucleus 

(fir' for the extended nucleus. 

r — r'\ 



The nuclear charge density pnuci(^) is defined by the nuclear model. In this paper we will use the Fermi 
model for the nuclear charge distribution. 

Here and in what follows we consider only the electric part of the classical electromagnetic interac- 
tion between the electron and the moving nuclei neglecting the magnetic interaction (e/c A(r)), which 
is small for low-energy collisions. 

The two-center expansion of the stationary wave function ^ri(^) and the time-dependent wave func- 
tion t) can be written in the form 

(5) 

*(^,^) = E ECaa(t)<^a,a(r-i?„(t)), 

a=A,B a 

where index a = A,B enumerates the centers, index a enumerates basis functions at the given center, 
and '^a,a{f — Ra) is the central-field bispinor, centered at point a. The coefficients c^^ of the expansion 
([5]) for the stationary wave function ^/'„ (r) can be obtained from the generalized eigenvalue equation 

Y.H,kcl = enY.S,kcl, (6) 

k k 

where indexes j and k enumerate the basis functions of both centers, and the matrix elements of H and 
S are given by 

Hjk = {j \hr,\k) , Sjk = {j I k) . (7) 



The expansion coefficients Caa{t) of the time-dependent wave function \l/(r,t) can be obtained by 
solving the linear system of first-order differential equations 

dCk{t) 



^ E ^^^—^ = E iH,k - T,,) at)) . (8) 

k k 

The matrix elements of T are given by 

TJk = ^{J\^\k)=T*,^+^^^S,,. (9) 

Obviously the matrix T is non-Hermitian, if the overlapping matrix S depends on time. 

The functions cpa depend on time due to two reasons. First, the basis functions centered at the target 
and projectile nuclei move together with the nuclei. Second, the basis functions depend parametrically 
on the distance between the nuclei, since their radial parts are obtained from the radial equations, where 
for each center the potential of the other nucleus is included in the so-called monopole approximation 
(see section HTBI) . Therefore, the time derivative of the basis function can be divided into two parts 



di 



dR 



dip. 



OR 



rfF ir^) - ^"/c ■ {^j I V I <^fc) , (10) 



where Va = dRa/dt is the velocity of the ion a. 
2. Trajectories of nuclear motion 

In the ion-ion collisions the intemuclear distance vector R = Rb — Ra, the length R of vector 
R, the target velocity (va), and the projectile velocity (vb) are time dependent. This dependence 
is defined by the trajectories of the nuclear motion. In low-energy collisions the nuclear trajectories 
can be obtained by solving classical non-relativistic Newton's equations of motion. In the case of 
point charges this solution is well-known Rutherford hyperbola (see Fig. [T]), which can be given in the 



parametric representation by the equations n25\\ 



where ^ G {—oo, oo) 



R = a{e cosh ^ + 1) 

t = - {e sinh ^ + ^ 
vo 



(11) 



i2 



(12) 



FIG. 1: The hyperbohc Rutherford trajectory (b is the impact parameter, Rq is the minimal distance between 
target A and projectile B, vq is the initial projectile velocity). The coordinate system St = {x, y, z) is defined 
with respect to the moving target ion. 

f is the initial velocity of the projectile, h is the impact parameter, and is the reduced ion mass. 
In the coordinate system = {x,y,z), which is shown in Fig. [H the X and Z components of the 
internuclear distance vector R are given by 



R cos 9 



X = Rsine 



where 9 = 2 arctg 



v/i^(th(e/2) + l) 

(£+l)-(£-l)th(e/2)^ 



(13) 



The angle 9 is related to the scattering angle 6oo by 6^ = ir — 9{t = oo) 



3. Time-dependent matrix Dirac equation 



In this work the two-center basis set Lpj is not orthonormal. Let us consider the transformation of the 
basis set ipj to the orthonormal basis ip^ by a matrix L^^ 



"Pi 



Vk 



^kj Vk ■ 



(14) 



Then the positive-defined matrix S can be represented as the product of L+ and L: 



S = L+L, Sf,= {^f I y.^) = {L-'^SL-'). = 6,,, . 



(15) 
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If the matrix L is an upper-triangle matrix, then the decomposition (1151) is so-called Cholesky factoriza- 
tion [1671]. The expansion of the time-dependent wave function over the orthonormal basis is given 
by 

^(^,^)= E E^iaW<a(r--i2aW,t), (16) 



a=A,B a 



where = LC. 



The time-dependent Dirac equation in the basis ipj' can be written in the form 

dt ^ ^ ' 

where M = — T^, the Hermitian Hamiltonian matrix is 



(17) 



Ht, = {^i\H\^,) = iL-'^ HL- 



(18) 



and the matrix is defined by 



Tl^^ = {^f\f\^f) = {L''^TL'%+tiL 



dL 



dt 



T-iL 



dL 
It 



(19) 



It should be noted that matrix is Hermitian, in contrast to the matrix T defined by Eq. Therefore, 
the matrix M is also Hermitian. 

The time-dependent matrix equation (fTTl) can be considered as a linear system of the first-order 
differential equations at the range of time t E (—00, 00). We assume that at the initial moment of time 
(t — > —00) the electron is localized on the target in the Is state and the projectile is the bare nucleus. 
Then, the wave function \l/(r, t) at t — )■ —00 is given by 



^(r,t) 



t— 00 



ipisir). 



(20) 



If the Dirac Is-target wave function i^isir) is included in the basis set, the initial conditions for the 
expansion coefficients can be written as 



t— >— 00 



t— >— 00 



(21) 



Equation (fTTl) is solved numerically, using approximate evolution operator 



C^(t + At) = e-'^'' C^{t) + 0{AH) 



(22) 



9 



where Hermitian matrix M is chosen as 

M = M(t + At/2). (23) 

Since the approximate evolution operator U{t) = exp {—iMt) is unitary, the time-dependent wave 
function conserves the norm at each time step 



{^{r,t)\^{r,t)) = Y.\Cf{t)\' = l. (24) 

j 

The matrix e"**^* is calculated at each time step using the eigen decomposition of matrix M 

M = VAV+, (25) 
where A is a diagonal matrix and columns of matrix V are the eigenvectors of M. Then one obtains 

^-iMt ^ y^-iAty+_ ^26) 

The time grid points U are chosen as ti = a/vo (e sinh + ^j), where the parameter ^ runs a uniform 
grid. The grid points Ri can be obtained using equation (fTT)) . 

B. Basis functions 

In our approach the basis set contains Dirac and Dirac-Sturm orbitals. The Dirac-Sturm orbitals 
can be considered as pseudo-states, which should be included in the basis to take into account the 
contribution of the positive- and negative-energy Dirac continuum. Both types of basis functions ifaa 
are the central field Dirac bispinors centered at the position Ra (a = A, B) 



r 



V 



(27) 



r 

where PnK{r) and QnK.{r) are large and small radial components, respectively, and k = {—lY'^^'^^^'^ij + 
1/2) is the relativistic angular quantum number. The large and small radial orbital components are 
obtained by solving numerically the Dirac or Dirac-Sturm equations in the central field potential V{r). 
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The Dirac equation is given by 



'-riK QuK 



(28) 



The radial components of the Dirac-Sturm orbitals which we denote by PnK{r) and Qn^ir) are the 
solutions of the Dirac-Sturm generalized eigenvalue equation 



+^]Qn.ir) + (v{r) + c' 



dr 



PnJr) 



A„.lV.(r)Q„,(r) 



(29) 



Here A„k can be considered as the eigenvalue of the Dirac-Sturm operator and WK{r) is a constant sign 
weight function. The energy is fixed in the Dirac-Sturm equation. If W{r) — i- at r — )■ cxo, all 
Sturmian functions have the same asymptotic at r — )■ 00. It is clear that for A„k = the Sturmian 
function coincides with the reference Dirac orbital which has the radial parts Pnonir) and QnoKi^). 
The widely known choice of the weight function is W{r) = 1/r, which leads to the well known 
'charge quantization' Z*^ = Z + A„k- The main advantage of this choice for the Coulomb potential 
V{r) = —Z/r is that the Coulomb-Sturmian orbitals can be given in an analytical form. This is not the 
case, however, for the non-Coulomb potential V{r). In the relativistic case the choice W{r) = 1/r is 
not very successful, because of the incorrect behavior of the Coulomb-Sturmian orbitals at r — )• 0. For 



this reason the standard form of the equation has to be modified [|5o[ 
In our calculations we use the following weight function 



68 



WJr) 



1 - exp(-(a«,r)2) 



(30) 



In contrast to 1/r, this weight function is regular at origin. It is well-known that the Sturmian operator 
is Hermitian and does not contain continuum spectra, in contrast to the Dirac operator. Therefore, 
the set of the Sturmian eigenfunctions forms the discrete and complete basis set of one-electron wave 
functions. 

The central-field potential V{r) in equations (|28] ) and (|29l) is arbitrary, and, therefore, it can be cho- 
sen to provide most appropriate Dirac and Dirac-Sturm basis orbitals. At short intemuclear distances 
the wave function of the electron experiences the strong Coulomb field of both nuclei. To take into 
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account this effect we have included the Coulomb potential of the second ion in the total one center po- 
tential V{r) in so-called monopole approximation. For instance, the total central-field potential V'^{r) 
of the center A is given by 

V^{r) = Kuci(r) + Kfon(0 , (31) 

where Vj^^i{r) is the Coulomb potential of the nucleus A and V^^^{r) is the spherically- symmetric part 
of the reexpansion of the potential V^uci(''' ~ ^b) on the center A 

Vlnir) = ^ Jd^A Kuci(^ - Rb) . (32) 

For the point nucleus the potential V^^^{r) is given by 

[ r>R 
VZ.{r) = J (33) 

C. Two-center integrals 

The matrix elements of H and S (Eq. |7]) are easily reduced to radial integrals f/O^, which are 
calculated by numerical integration in the radial semi-logarithmic grid [711 . 

Modified Lowdin reexpansion procedure 

Two-center matrix elements are calculated using a symmetrical reexpansion procedure, proposed in 



Refs. [156, 



57l1 . The reexpansion procedure is based on the technique developed by Lowdin [|55[1 . We 
assume that in the local coordinate frame the ^-axis is directed along the internuclear axis A-B (see 
Fig. O. The following geometrical relations take place 



va = r - Ra, vb = r - Rb , R = Rb - R/ 



2 I p2 „2 „2 d2 „2 (34) 

n ■ A ' ^ B n 

COS OA = :Z7^ , COSUb 



2RrA 2RrB 
Let the indexes a and b enumerate basis functions centered at the points A and B, respectively. The 
standard Lowdin reexpansion of the nonrelativistic central-field function Fb{rB) centered at the point B 
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FIG. 2: Integration regions Va and Ve- 



in terms of the spherical harmonics 11m (^"a) centered at the point A can be written in the form [155 . 



72d 



Fb{rB) = ^^^^ Yi^^^^iOn, = — J2 o^irnMb, lb\rA) Yim^iOA, ^) , (35) 
tb Ta 

where aim^{fb, /fe|r^) is so-called Lowdin a-function defined by 

airnAhMrA) = ^^^^^^ J Mr) P^' '^f^' '1 Pj--' j dr . (36) 

Here P^'™' is the standard associated Legendre polynomial and Kim is the normalization constant 



Im 



2l + l{l-\m\)\ 

\ 2 {i + \m\y: ^ ^ 

Similarly, the function Fa(r^) centered at point A can be expanded in spherical harmonics YimiOs, ^) 
centered at the point B. 

When the logarithmic or semi-logarithmic grid is used the radial grid step increases with increasing 
radius ta- Therefore, the Lowdin reexpansion procedure becomes unstable and poorly convergent for 
the values of radius ta in the region near ta = R, especially for the oscillating and strongly localized 
atomic-like wave functions. In addition, the Lowdin procedure is not symmetric with respect to the 
centers A and B. 

To improve the convergence we modified the standard Lowdin reexpansion procedure by dividing 
the range of the integration into two regions Va and Vb as shown in Fig. [2l The region Va contains ion 
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A and the region Vb contains ion B. The dividing of the integration area into two parts can be done, for 
example, by a plane passing through the center of the segment (AB). We apply the reexpansion proce- 
dure only to the "tails" of the wave functions occurring in a given region. To describe this procedure 
we introduce the step-wise functions 0yi(r) and ©^(r) by 



1 reVA 
veVb 



r eVa 

1 teVb, 



(38) 



and rewrite the product of the functions centered at the different points in the following way 

Fair a) ■ F,{rB) = F„(r^) ■ {F,{rB) G^W) + e^W) ■ F,{rB) . 



(39) 



The reexpansion of the function tail Fb(rB) ©aI^') centered at B onto center A has the form 

FtivB) ©aW 

where 



^^Yi^mS^B,^)eA{r) = — J2^imMbJb\rA) ■Yi„,^{9a,v), (40) 
tb Ta I 



aimSfb,h\rA) 



R 



\rA+R\ 



Mr) Pi 



\mb\ ( ^A 

b 



2rR 



I 2rAR 



dr (41) 



and = max{rA, — 



In the relativistic case the spin-angular part Xkui of the large and small components of the central- 
field wave function is the Pauli spinor [73] 



(42) 



m,ms ' ^ ' 

where C-j^ i are the Clebsch-Gordan coefficients il4\ and (a) is a spin function. 

The symmetric reexpansion of the relativistic wave function "tails" onto centers A and B can be 
written in the form 



I — XK,^,^rB) \ 

rB 

\ Tb ' 



&A{r) = y: 



X^^,t(fA) 



TA 



■rA 



\ 



and 



rA 

eB(r-) = E(-i)'"-^ 



X-K^ib{rA) 



(43) 



TA 
Qa{rA) 



I « X-Kaila\TA) ] 

\ TA / 



TB 



TB 



XkiiA^b) 

X-^^ilaifs) 



(44) 
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The p- and g-functions, which are the relativistic analogs of the modified Lowdin a-functions, are 
defined by 

(45) 

mi,, ma 

where 1 = 1 — sign(K). Functions aim,,{PbJb\rA) and a-i^^{QbJb\rA) are defined by equation (|4TI) . 
where the function fb{r) has to be replaced by the functions Pb{r) and Qb{r), respectively. Functions 
PK^lai^VB) and qKi_ia{^VB) are defined similarly to equation (l45l) . where indices A and 6 should be 
replaced by B and a, respectively. 



Two-center overlap integrals 



Let us consider two-center overlap integrals S^^^ . Here and below symbol (0) means that the integral 
is considered in the local coordinate frame, where the z-axis is directed along intemuclear axis A-B. 
The integral 5*^°^ can be divided into two parts 



S. 



= (a I b) = {a I b)A + {a \ b) 



B 



(46) 



where the notations <>a and <>b mean the integration over the regions Va and Vb, respectively (see 
Fig. [21). Using the reexpansions of the large and small components onto the center A (in the region Va) 
and onto the center B (in the region Vb) we obtain 



(a I b)A 
(a I b)B 



Pa{r)-p,^^Sb\r) + Qa{r)-q,^^Sb\r) 



(47) 



Matrix elements of the nuclear attraction potentials, V^^i^ta) and Kmci(^B)' and of the mass operator 
/3mc^are calculated similarly to the overlap integral. 
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Two-center gradient matrix elements 



As in case of the overlap integral, the region of integration for the gradient matrix element (g) 
is divided into two parts, 



G^aho) = (a I V, I 6) = (a I V, I 6)^ + (a | V, | 6) 



B ■ 



(48) 



Here index q = 1, 0, —1 enumerates covariant spherical coordinates. Using the Gauss theorem nl5n for 
the integration over region A, we obtain 



(49) 



where {a \ b)s is the surface integral over the region S (see Fig. O. The volume integrals over regions 
A and B are given by 



(50) 



oo 

(H V, I a)A = E a^'Uf^b 

K Q 

oo 

{a\Vq\h)B = 5l'(-l)'""V'^0>a,j6/Uf>) j dr^^^^{a\r)b^^^^Pb{r)+q^^^{a\r)b^^^^^Mr) 

K g 

where the prime at the sum symbol indicates that the summation is restricted to odd values of la + 1 and 
lb + I, and the operator -Dk,k' is defined by 



d + 1) - k{k + 1] 

k' ~ ~; — I 

dr 2r 



(51) 



The coefficients g^'^H^i-, j V') the relativistic analogs of the Gaunt coefficients TI6\ 



,/ 2j + l 2j' + l 



(52) 



The relativistic Gaunt coefficient is non-zero only if / + /' + A; is even. 
The surface integral is given by 

oo 

(a \b)s = S. -\y C'^^\ C''^\ I dr- PJr) PAr) Ui uj— 



R/2 



r 1 /? 
X y Cf / rfr - Qa(r) Qdr) U-, -, J —) 

m,ms 



(53) 



where 



Ui^,i,M = {-l)'^-^^{2U + l){2h+l)Ki^\rn\Ki,\„,\PC\x)P, 



'b 



(54) 
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Two-center (a. ■ p) matrix elements 

Two-center (a ■ p) matrix elements A^^^ can be divided into three parts, similarly to the gradient 
matrix elements G^^^{q) (l49l) . 

A^^^ = {a\oL-p \ h) = {h\oL- p\a)A + (a | o: • p | 6)b + - (a | ckq | 6)5 , (55) 
where the volume integrals (6 | a ■ p | a)^ and (6 | a ■ j? | a)^ are given by 



00 



(56) 

The last term in equation (l55l) is the surface integral, which is given by 

1 7 1 

-{a I «o I 6). = E -^^7^:4..^/::^. y dr-Q^ir)P,{r)Uj^,^^{R/2r) 

00 

+ E / dr- P^ir)Q,ir)U,j^^{R/2r). 



(57) 

rfr^P,(r)g,(r)[/,j^„(i?/2r). 
Transformation of the two-center matrix elements to the laboratory frame 

As indicated above, the laboratory frame S is defined to be at rest with respect to the initial target 
position. Then, the two-center matrix elements calculated in the local coordinate frame S' = (x', y', z') 
(see Fig.[T]) have to be transformed to the laboratory frame S = (x, y, z). The corresponding two-center 
integrals can be obtained from S^^^^ , and G^^f^ by rotating the coordinate system around the y- 
axis for angle —9 (Fig.[T]). For the overlap integrals Sab and the two-center Dirac-Hamiltonian matrix 
elements Hat in the laboratory frame S we obtain 

Sa,b = SnaKatJ.a;ni,Ki,tii, ~ E '^MaM^-^) dfi,,fj,{6) Sli^i^^^. ji^i^^^ , 

^ J* (58) 
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where di'^i^iO) are real Wigner's D-functions uM- The transformation of the gradient matrix elements 
G^^^ (q) to the laboratory frame S is given by 

D. Charge-transfer probabilities and cross sections 

1. Transition amplitudes 

Transition amplitude for electron capture to an ion state an is given by 

T„„(t) = {tfJanir, t) I ^(r, t)) , t ^ oo . (60) 

As previously mentioned, here index a = A,B enumerates different centers (target and projectile ions) 
and ipanii", t) are the wave functions of the free-moving ion a. After the collision (t oo) the wave 
functions ^/'^^(r, t) of the free-moving ion a are given by 

ij^nir, t) = e-^^""* s^{r) < Jr - R^) , (61) 

where i'ani''^) the stationary Dirac wave functions of ion a at the rest and 50,(7") is the translation 
factor. For the low-energy collisions the translation factor and the energy Fan of the moving ion can be 
taken in the nonrelativistic approximation 

So(r) = exp{iVa ■ r) , E^n = £an + vl/2 . (62) 

Generally, the translation factor is introduced in the basis functions to improve the convergence of the 
time-dependent wave function expansion. We did not include the translation factor Sair) in the basis 
functions because of the computational complexity of the two-center integral calculations. However, 
at the limit t — )• 00 we can reexpand the moving orbitals Saffan over the basis (fan in the following 
way. At the limit t — )■ 00 the basis functions of the different centers do not overlap. The basis set Lpj is 
orthonormal and 
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The expansion (|63|) is not exact, since the finite basis ^pan is incomplete and K is non-unitary matrix. 
In particular, 

I Kan', an? < 1 • (64) 



We can renormalize matrix K and rewrite expansion (1631) in the form 

e-'^^-^^nW - Y.Kan',an^an'{r), K = NK, N = {K K+)-'/\ (65) 

n' 

Here matrix plays a role of the normalization factor. The renormalized matrix K is unitary. 

Now we can obtain the expansion of the time-dependent wave function \E'(r, t) over the basis func- 
tions with translation factor (t — )■ oo), 

*(r, t)=Y. cut) ^a,n{r) - E Canit) e*""''' ^anir) , (66) 

a,n a,n 

where 

Can{t) = E Kan,an' Can' if) ■ (67) 
n' 

It should be noted that the set of coefficients Can{t) is normalized to unity 

E|C'an(t)r=E|Can(t)|' = l, t OO . (68) 

Therefore, for the transition amplitude Tan{t) we obtain 

Tan{t) = {i^anir, t) | ^(r, t)) = E Can'{t) c'^^' | ^an') , t ^ OO . (69) 

n' 

The stationary Dirac wave functions including wave functions of the positive and negative energy 
spectra, form a complete basis set. Therefore, 

EI^-I' = E l^-WI' = i- (70) 

Q,n Q,n 

2. Transition probabilities 

Transition probabilities Wan{t) are defined by 

Wanit) = \Tan{t)?. (71) 
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The probability Wanit), defined by equations (TtTI) and (l69l) . has an oscillatory behavior at t — )■ oo, 
because the basis functions Lpa,n are not the solutions of the hydrogen-like Dirac equation and the basis 
set is truncated. We can remove the oscillatory component of the probability Wan{t) for the large time 
(t —7- oo) in the same way, as it was done in Refs. tA\m- 

At the large time (t oo) the coefficients Cj'(t) coincide with the coefficients Cj{t). Therefore, 
the coefficients Cj (t) are the solutions of equation (flTI) 

' ^ = ^^^^ ^^^^ • ^'^^^ 
Then for the coefficients Cj (t) we obtain the equation 

I '^^^ = KMit) Cit) . (73) 
at 

Using the diagonalization procedure for the Hermitian matrix K M , we can decompose 

KMK^ = VVlV^, (74) 

where ^7 is a diagonal matrix with eigenvalues Vtkk = and V is a unitary matrix. 
We introduce new coefficients B{t) by 

Bu{t) = (e^^*r+C(t))^ = e-'=*Ev;.;C,(t), J2\Mt)\' = l. (75) 

j k 

These coefficients have well defined limits at t — )■ oo. The amplitudes TQ,„(t), defined by equation (|69|) , 
in terms of the coefficients B{t) are given by 

Tj{t) = e^(^^--'=)* V^fc B, I ^i) . (76) 

k,l 

Then for the probabilities Wj (t) we obtain 

Wj{t) = \Tj{t)\^ = ^ \Bk{t)r 1^ Vik{^'^ + (oscillating term) . (77) 



|2 

k ' ' I 

Removing the oscillating term (S'], we can introduce probabilities Wj{t) defined as 

|2 



w;{t) = j2\BM \ ^i) . (78) 



20 



Since hydrogen-like Dirac wave functions i})^ ^ of each center a (including the positive and negative 
Dirac continuum spectra) form a complete basis set, we get 

E mt)? = E \m\' E \v^i\' = E \m\' = i • (79) 

j&a k I k 

The coefficients -Bfc(t) and the matrix elements have well-defined limits for t oo. Therefore, there 
exists the limit 

P^,n = lim W^t) . (80) 

t— ^00 

The direct (Pd), charge transfer (Pet), and ionization (Pion) probabilities are given by 

Pd = E } -Pet = E ^B,n , Pion = ^ ~ E "^Cf" = ^ ~ Pd — Pet • (81) 



a,n 



where the prime at the sum symbol indicates that the summation runs over the discrete bound states of 
the ion a. 

The cross sections for the charge transfer (dct) and ionization ((Xion) processes are then calculated as 
usual by integrating the probabilities over the impact parameter b: 

oo oo 
(Tet = 27f J db bP,t{b) , aion = 2n J db bPionib) • (82) 



3. Z-scaling 

It is well-known that in the nonrelativistic theory the scale transformation r' = Z r and R' = Z R 
allows one to transform the wave functions V^(r) and the energies £ of a homonuclear one-electron 
quasi-molecule with a point nuclear charge Z > 1 to the wave functions ip'ir) = Z^^/'^il){Zr) and 
energies e' = e/Z'^ of the molecule. The same scale transformation can be considered in the 
nonrelativistic homonuclear collisions. 

The time-dependent Schrodinger equation describing A^^^^^^-A^^ collision is given by 

(1 Z Z \ d 

2 \r-RAit)\ \r-RBit)\J ^ ' ^ dt ^ ' ^ ^ ' 

If in Eq. ([83]) we set r' = Z r, t' = Z'^t, and = Z Ra{t), we obtain the time-dependent 



Schrodinger equation for the H+-H collision y. 



121 
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It should be noted that the scaling R! = ZR(t) is satisfied exactly for a straight-line trajectory. In this 
case the impact parameter b, the velocity Va, and projectile energy E are transformed by 

b' = Zb, v'^ = Va/Z, E' = E/Z\ (85) 

It follows that the probability Pq,,„ (Z, b, E) and the cross section aa,n{Z, E) for a process in a symmetric 
system with nuclei of charge Z can be obtained from the probability Pq „(1, 6', E') and the cross section 
(Tq,,„(1, E') for the same process in the H(ls)-H+ system by the relations 

P„,„(Z, 6, E) = P,,„(l, b Z, E/Z^) , E) = ^ E/Z^) . (86) 

It should be noted that this scaling law is not valid for the relativistic collisions. 



III. RESULTS 



A. The choice of the basis 



In our relativistic calculations, we used two Dirac-Sturm bases of different size. Both bases in- 
clude functions of the positive-energy Dirac spectrum and Sturm orbitals corresponding to the negative- 
energy Dirac spectrum. It should be noted that the constructed bases satisfy the dual kinetic balance 
conditions itzzl] and do not contain so-called "spurious" states [78]. 

The positive-energy functions of the first basis (Basis 1) on each center in the standard nonrelativistic 
notations of atomic shells are given by: ls-3s, 2p,3p, 3d, 4s-6s, 4p-6p, 4d-6d, 4f,5f . Here the overline 
symbol (nl) is used to indicate the Dirac-Sturm (pseudo-state) basis functions. The total number of the 
positive-energy orbitals of both centers is 220 and the total size of Basis 1, including the negative-energy 
orbitals, is 440. Basis 1 is used for both stationary and time-dependent wave functions. 

The basis size can be increased in the calculations of the stationary states of quasi-molecules. 
Positive-energy functions of the second basis (Basis 2) is constructed from 26 atomic shells: Is, 2s-8s, 
2p- 8p, 3d- 8d 4f- 6f 5g, 6g. In Basis 2, the total number of orbitals of both ions, including the 
negative-energy spectrum, is equal to 784. This basis is used only in the calculations of the stationary 
states of quasi-molecules. 
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B. Stationary ground states of some homonuclear quasi-molecules 



1. Energies of the ground state of some homonuclear quasi-molecules 



In Table HI we present the results of our relativistic calculations of the lag state energy of the H^, 
Th2'^^^, and U2^^^ quasi-molecules for so-called chemical distance R = 2/Z a.u.. Since the quasi- 
molecule Thg^^^ was considered as a reference system for testing relativistic effects, it was calculated 



in a number of papers using high-precision large-scale methods (see, e.g., Refs. 



As one 



TABLE I: Relativistic energies (a.u.) of the lag state of quasi-molecules for the point-charge nuclei and R = 2/Z 
a.u.. 





H+ {Z = 1) 


Th^^9+ {Z = 90) 








Rel. error 




Rel. error 




Basis 1 
Basis 2 
Others 


-1.1026248 
-1.1026405 
-1.1026416'" 


1.5 • 10-^ 
1.0 • 10-6 


-9504.573 
-9504.732 
-9504.756'' 


1.9 • 10"^ 
2.5 • 10-6 


-9965.190 
-9965.307 



Ref . [58J 
^ Ref. [59J 



can seen from TableU there is a good agreement of our data with very accurate values obtained in Refs. 

The relative precision of our results for the quasi-molecules and Th2'^^^ is increased by an 
order of magnitude when Basis 1 is replaced by Basis 2. 

In Fig. [3] we display the energy of the lag state of the U2^^^ quasi-molecule as a function of the 
internuclear distance i? on a logarithmic scale. In this figure the solid line indicates the energy E(R) 
calculated using two-center Dirac-Sturm Basis 2. The dashed line represent the results of the one-center 
calculations in the monopole approximation. As one can see from Fig. [3l in the two-center basis the 
lag electron "dives" into the negative-energy Dirac continuum at a critical distance Rc = 34.7 fm. 
The critical distance obtained in our monopole approximation amounts to i?c = 25.5 fm that is too far 
from the exact value. It should be noted that in our one-center monopole approximation the basis was 
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91+ 92+ 

Energies E (in relativistic units) U -U 



— Is (monopole approxmation) 
la (two-center basis) 



R,(ls)=25.5 fm 
R.(la)=34.7 fm 




R,(ls) R^(lo) 



R[fm] 

FIG. 3: The lag state energy of the U2 quasi-molecule as a function of the intemuclear distance i? on a logarith- 
mic scale. 



centered at the position of the nuc 



eus (A or B) but not at the center of the internuclear interval, as it 



was done in the papers |141|,|62|,|82D. The monopole approximation used in Ref. HlW is more suitable 



for the short-distance regime and gives the value of the critical distance equal to Rc = 35 fm. 



2. Critical distance 



In the Table HI] we present our results of the two-center relativistic calculations (Basis 2) of the 

(2Z 

critical distance -Rc for homonuclear one-electron quasi-molecules A2 and compare them with the 
corresponding values obtained by other authors. There exists a discrepancy of about 5-10% between 



the critical distance data for the point nuclei [1271 . 141 
in a very good agreement with the results of Ref. 



m. 



60, 



61 



83L 



8411 . Our results for this case are 



In our calculations for extended nuclei, the 
finite nuclear size was taken into account using the Fermi model of the nuclear charge distribution (for 
details, see, e.g., Ref. [i85J). The root-mean-square nuclear charge radii (-R^)^^^ were taken from Refs. 
186] (for Z=88), ifl (for Z=90), [88] (for Z=92), and [89] (for Z=9A, 96, 98). The number of works 



where the finite nuclear size effect was taken into account is much less than for t 



le point nucleus case. 



We can systematically compare our results only with the data obtained in Ref. []63]l . The discrepancy 
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TABLE II: Critical distances (fm) for homonuclear one-electron quasi-molecules A2 



(2Z-1)H 





Point nucleus 


Extended nucleus 


z 


This work 


Others 




This work 


Others 


88 


24.27 


24.24'' 


5.5705 


19.91 


19.4'= 


90 


30.96 


30.96*^ 


5.7210 


27.05 


26.5'= 


92 


38.43 


38.42" 


5.8569 


34.72 


34.3'= 






36.8^ 






34.7'^ 


94 


46.58 


46.57'' 


5.794 


43.16 


42.6'= 


96 


55.38 


55.37*^ 


5.816 


52.09 


51.6^ 


98 


64.79 


64.79'" 


5.844 


61.63 


61.0^ 












61.1'^ 



'' Ref. [60]. Ref. [27] > Ref. [63], Ref. [62] 



between our data and those from Ref. 16311 is considerably larger for the extended nuclei than for the 
point nuclei. A possible reason of that could consist in a rather crude estimate of the nuclear size effect 



in Ref. 16311 . It should also be noted that in the work 



63D other values of the nuclear radii, namely 
3/5 ■ 1.2 ■ (2.6 • Zy/^ were used. Our calculations showed, however, that the usage of the 
nuclear radii from Ref. ll63n changes the values of not more than by 0.02 fm. 



C. Charge-transfer probabilities and cross sections 

1. H(ls)-H^ collisions 

Fig. m shows the charge-transfer probabilities Pct{b) for the H(ls)-H+ collision as functions of 
the impact parameter h for the projectile energies of 2 keV and 5 keV. The results of our relativistic 
calculations for 2 keV (solid line) and 5 keV (dashed line) are found to be very close to nonrelativistic 
calculations based on the two-center atomic-orbital (AO) expansion [15i|90|l. This is not surprising, since 
the relativistic effects are negligible for the H(ls)-H+ collision. Our calculations were performed for 
the straight-line trajectory of the projectile that corresponds to the full screening of the target nuclear 
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3 4 

b (a.u.) 



FIG. 4: Charge-transfer probabilities Pch(b) for the H(ls)-H+ coUision as functions of the impact parameter b. 
Our resuhs for the projectile energies 2 keV (solid line) and 5 keV (dashed line) are compared to the related 



results from the AO-expansion calculations 



M (symbols "+" and "x"). 



charge by the Is-electron. 

As was demonstrated in Ref. [5*], the two-center AO expansion data are in a very good agreement 



with results obtained by a direct numerical solution of the nonrelativistic Schrodinger equation 



In the review |15|], the two-center AO expansion data 119011 are also compared with the results 1191 . 



|14]- 



obtained by the expansion over the "nonmoving" Hylleraas functions. It should be noted that the 



9011 . [|14n and to our results in trend but 



Hylleraas expansion data are similar to the results of Refs. 
differ by the phase and the magnitude. 

In Table Hill we present the total charge-transfer cross sections for the H(ls)-H+ collision in a 
wide range of the projectile energy (from 0.5 keV to 100 keV) and compare them with nonrelativis- 
tic large-scale calculations of the recent paper [|l3ll . which can be considered as an extension of the 
pioneering works [18i[93|], where the analytical Sturmian basis set expansion was used. We also give 
the cross section values, deduced from the experimental results n9m . The relative uncertainties of the 
recommended and interpolated experimental data are about 5%-10%. As one can see from the table Hill 
our results are in a good agreement with the theoretical data of Ref. /l3^ and with the experimental 
data. 
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TABLE III: Charge transfer cross section for the H(ls)-H"'" coUision, in units of 10 cm^. 



Projectile energy 








E (keV) 


This work 


Winter [13J 


Exp" 


0.5 


199.6 






0.7 


186.9 






1.0 


172.4 


173.0 


171 


2.0 


144.9 




144 


4.0 


117.5 


118.1 


115* 


5.0 


107.8 




110 


10.0 


81.3 




77.5 


15.0 


63.5 


67.41 


55.6* 


20.0 


48.9 




44.5 


25.0 


36.2 


39.45 


35.3* 


30.0 


26.6 




27.6* 


40.0 


15.3 




16.5* 


50.0 


9.1 


10.04 


9.9 


60.0 


5.6 




5.9* 


70.0 


3.5 




3.6* 


80.0 


2.3 




2.3* 


100.0 


1.1 


1.11 


1.1 



" Recommended values [941 deduced from the experimental data. 
* Interpolated values obtained using an anahtyc fitting function [941 . 



The ionization cross sections, computed in this work using equations (|8T1) and (|82| ). are displayed in 
Fig. [51 Our results are in a good agreement with the experimental data in the range of the proton energy 
from 20 keV to 80 keV. At the energies less than 15 keV we observe a significant relative deviation of 
our results from the experimental data. Thus is probably due to the fact that the absolute uncertainty of 
our data is approximately the same in the whole region of the energies ( about (l-3)-10^^^ cm^), while 
at the low energies the ionization cross section tends to zero. This leads to a large relative error in the 
low energy region. 
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FIG. 5: Ionization cross section o"ion(E) for the H(ls)-H"'" collision as a function of the projectile energy E. The 
solid line is obtained by the interpolation of our results indicated by circles, the dashed line is obtained by the 



interpolation of data from Ref. 11311 . and the triangles indicate the experimental data from Refs. ||95 



m. 



n 

In contrast to our results, the theoretical data of Ref. [13], which are shown in Fig. |5Jby squares, 
are in a good agreement with the experimental data in the low energy region (less than 25 keV) and 
significantly differ (at least by 25%) from the experimental data for the energies larger than 40 keV. The 
reason of this discrepancy is unclear to us. 



2. Ne^^(ls)-Ne^^+ collisions 

To study the role of the relativistic effects in the homonuclear collisions and to test our approach we 
calculated the charge transfer cross sections for the Ne^+(ls)-Ne^°+ collisions with the standard value 
of the speed of light (c = 137.036 a.u.) and in the nonrelativistic limit (c — )■ oo) by multiplying the 
standard value of the speed of light by the factor 1000. The obtained values are presented in Table HVl 
It should be noted that the projectile energy values are divided by Z"^ {Z is the nuclear charge) and 
the values of the ionization cross section cTjon are multiplied by the factor Z"^. This was done in order 
to compare the Ne^+(ls)-Ne^°+ cross section data with the H(ls)-H"'' results in accordance with the 
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TABLE IV: Charge transfer cross section act{E) (10 cm^) as a function of the projectile energy E for the 
Ne9+(ls)-Neio+ (Z=10) and H(ls)-H+ collisions. 





Ne9+(l£ 


)-Neio+ 


H(ls)-H+ 


Ne9+(ls)-Neio+ 


E/Z^ 




actiE)-Z^ 


(Jct{E) 


actiE)-Z^ 


(keV/u) 


Rel.'* 


Nonrel.^ 




Born approximation'^ 


1.0 


171.6 


172.2 


172.4 


188.4 


2.0 


144.3 


144.8 


144.9 


150.7 


4.0 


117.1 


117.5 


117.5 


114.8 


5.0 


107.3 


107.7 


107.8 


107.3 


10.0 


80.8 


81.3 


81.3 


76.2 


1 ^ n 








J / .u 


20.0 


48.5 


48.9 


48.9 


48.2 


25.0 


35.9 


36.2 


36.2 


38.1 


30.0 


26.4 


26.7 


26.6 


30.1 


40.0 


15.1 


15.3 


15.3 


19.9 


50.0 


9.0 


9.1 


9.1 


13.7 


60.0 


5.6 


5.6 


5.6 


9.1 


70.0 


3.5 


3.5 


3.5 


5.4 


80.0 


2.3 


2.3 


2.3 


3.6 


100.0 


1.1 


1.1 


1.1 


2.0 



Relativistic calculations. 
^ Nonrelativistic limit (c — )• oo). 
^ Born approximation [971 . 



scaling low (|86l) . As one can see from the table, the relativistic effects, which decrease the values of 
the charge-transfer cross section, are rather small and can be estimated as 0.5-0.8 (aZ)'^, where a is the 
fine-structure constant. In Table HVl we also compare our scaled nonrelativistic Ne^+(ls)-Ne^°+ data 
with the H(ls)-H+ results. It should be noted that our calculations for the Ne^+(ls)-Ne^°+ collision 
were performed for the Rutherford trajectory (see section lll A2I) . This is probably the reason of a very 
small discrepancy between the data presented in the 3-rd and 4-th columns of Table HVl 
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b{fm) 



FIG. 6: The charge-transfer probabiUty Pch{b) for the Xe^^+(ls)-Xe^^+ coUision as a function of the impact 
parameter b. The sohd line interpolates the relativistic values (squares) and the dashed line corresponds to the 
nonrelativistic limit. In both cases, the projectile energy is £'=5.9 MeV/u. 



It is also of interest to compare our results with the calculations performed within the plane wave 
Bom (PWB) approximation. The results of such a calculation for the Ne^+(ls)-Ne^°+ collision [|97D 
are presented in the 5-th column of Table HVl The details of the modified PWB method can be found 
in Ref. [19 811 . It is seen from the table, that the PWB data are in a reasonable agreement with the more 
elaborated calculation. 



3. Xe^^+(ls)-Xe^^+ collisions 

The relativistic effect for the Xe^^'^(ls)-Xe^^"*" collisions is considerably larger than for the 
Ne^+(ls)-Ne^°+ collisions. The computed relativistic (solid line) and nonrelativistic (dashed line) 
charge transfer probabilities Pct(&) as functions of the impact parameter b for the projectile energy of 
5.9 MeV/u are displayed in Fig. [6l The oscillatory behavior of both curves is the same but the nonrela- 
tivistic curve is shifted toward higher energies. 

In Table rvl we present the relativistic and non-relativistic values of the charge-transfer cross section 
for the Xe^^+(ls)-Xe^^+ collision scaled to Z = 1. One can see from the table, the relativistic effect 



30 



increases from 10% to 40% with increasing the projectile energy. 

TABLE V: The charge transfer cross section (Jct{E) (10^^^ cm^) as a function of the projectile energy E for the 
Xe^+(ls)-Xe^°+ and H(ls)-H+ colhsions. 





Xe53+(i5).xe54+ 


(Z=54) 


H(ls)-H+ (Z=l) 


E/Z^ 


E 


ME) ■ Z^ 


(Jct{E) ■ Z2 


^ct(^) 


(keV/u) 


(MeV/u) 


Rel.'^ 


Nonrel.'' 




1.23457 


3.6 


148.3 


163.3 


165.0 


2.02332 


5.9 


129.4 


143.0 


144.9 


3.42936 


10.0 


109.1 


123.8 


124.8 


34.2936 


100.0 


13.3 


20.6 


20.7 



" Relativistic calculations. 

^ Nonrelativistic limit (c — >• oo). 



4. U^^+(ls)-U^'^+ collisions 

The calculations of the charge-transfer probabilities and cross sections for the U^^+(ls)-U^^+ colli- 
sions were performed for the extended nuclei. The Fermi model of the nuclear charge distribution with 



-Rnuci = 5.8569 fm was used [18 811 



The computed relativistic (squares) and non-relativistic (circles) values of the charge-transfer prob- 
abilities Pct(^) and the interpolating curves are displayed in Fig. Ul It is seen from the figure, the 
nonrelativistic and relativistic probabilities significantly differ. The nonrelativistic curve (dashed line) 
is shifted toward higher energies compared to the relativistic one (solid line). The same curves in the 
small impact parameter region are shown in Fig. [81 In this figure, the vertical dashed line indicates 
the critical impact parameter be = 27.5 fm. For b = be the lag ground state level of the U2^'^^ quasi- 
molecule reaches the negative-energy Dirac continuum. It should be noted that for the non- straight-line 
(Rutherford) trajectory the value of the critical impact parameter be is less than the critical distance 
Rc presented in Table HIl For values b smaller than be the lag level "dives" into the negative-energy 
continuum. 
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b(fm) 

FIG. 7: Charge-transfer probability Pch(fc) for the U^^+(ls)-U^2+ collision as a function of the impact parameter 
b. The solid line interpolates the relativistic values (squares) and the dashed line corresponds to the nonrelativistic 
Umit. In both cases, the projectile energy is E = 6.0 MeV/u. 
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FIG. 8: Charge transfer probability Pet (b) for the U^^^(ls) - U^^"'" collision as a function of the impact parameter 
b in the small b region. The value b = be corresponds to the diving of the lag level into the negative-energy 
Dirac continuum. The solid line interpolates the relativistic values while the dashed line corresponds to the 
nonrelativistic limit. 
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In Table rvTl we present the results of our relativistic (3-rd column) and non-relativistic (4-th column) 
calculations of the total charge-transfer cross section actiE), scaled to Z = 1, for the U^^+(ls) - U^^"*" 
collision at different values of the projectile energy E. The values of cXctiE) were obtained for the 

TABLE VI: Charge transfer cross section act{E) (10^^'' cm^) as a function of the projectile energy E for the 
U9i+(ls)-U92+ and H(ls)-H+ collisions. 







H(ls)-H+ 


Energy 


Energy 


act ■ 


act ■ Z^ 


^ct • Z^ 


act 


E/Z"^ (keV/u) 


E (MeV/u) 


Rel. 


Nonrel. 


Nonrel. str. line 




0.70889 


6.0 


135.3 


184.2 


185.0 


186.4 


0.76796 


6.5 


132.7 


181.3 


182.0 


183.1 


0.82703 


7.0 


130.3 


178.2 


179.1 


180.1 


1.18147 


10.0 


117.1 


165.8 


166.7 


167.6 



Rutherford trajectories of the target and projectile ions. One can see from the table, the relativistic 
effect amounts to about 30% of the non-relativistic value of a^- In the 5-th column of Table IVTl we also 
present our results obtained for the straight-line trajectory of the projectile ion (in this case the target ion 
is at rest). As one can see from the table, the difference between the results obtained for the straight-line 
trajectory and the Rutherford one is very small. The non-relativistic values of the charge-transfer cross 
section for the U^^+(ls)-i-U^^+ collision, scaled to Z = 1, are also compared with the cross section 
cTctiE) for the H(ls)-H+ collision, presented in the 6-th column of Table IVTl 

IV. CONCLUSION 

In this paper we presented a new method for the relativistic calculations of one-electron two-center 
quasi-molecular system in both stationary and time-dependent regimes. The method is suitable for a 
wide range of the intemuclear distances including the critical regime, when the ground state of the 
quasi-molecule can dive into the negative-energy Dirac continuum. Using this method we calculated 
the energies of the H2, Thg^^^ and Ul^^~^ quasi-molecules, the critical distances for some homonuclear 
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quasi-molecules A+(^^ (Z=88, 90, 92, 94, 96, 98), and the charge transfer probabilities, charge trans- 
fer and ionization cross sections for the H(ls)-H+, Ne9+(ls)-Ne^°+, Xe^^+(ls)-Xe^^+, and U^i+(ls)- 
U^^+ low-energy collisions. 

The results of our calculations of the charge transfer probabilities and cross sections for the H(ls)- 
collision are in a good agreement with experimental data and with theoretical results obtained by 
other authors. The influence of the relativistic effect on the charge transfer probabilities and cross 
sections for the Ne9+(ls)-Nei°+, Xe^^+(ls)-Xe^^+, and \J^^+(ls)-\J^^+ collisions is investigated. We 
demonstrated, that the relativistic and nonrelativistic charge-transfer probabilities as functions of the 
impact parameter have the same oscillatory behavior at low energies, but the relativistic curve shifted 
to lower energies compared to the nonrelativistic one. In the case of the U^^^(ls)-\J^'^~^ collision the 
relativistic effect reduces the values of the cross section by about 30%. 

In our further investigations we plan to study in more details the effect of diving the lag level of the 
Ua^^"*" quasi-molecule into the negative-energy Dirac spectrum and the influence of this effect on the 
values of the charge-transfer probability. With this goal, we are going to develop an approach which 
would allow us to compare the calculated probabilities with and without the diving of the ground state 
into the negative-energy continuum. We also plan to extend our method to collisions involving many- 
electron ions and neutral atoms. This will allow us to study the Is-ls charge transfer in low-energy 
heavy ion-ator n col lisions. Such experiments, that were successfully performed for low- and middle-Z 
collisions [1991 - 1 10211 . are presently under preparation for high-Z collisions at GSI. 
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